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Abstract 

Most models for asymmetric dark matter allow for dark matter self annihilation 
processes, which can wash out the asymmetry at temperatures near and below the 
dark matter mass. We study the coupled set of Boltzmann equations for the symmet- 
ric and antisymmetric dark matter number densities, and derive conditions applicable 
to a large class of models for the absence of a significant wash-out of an asymmetry. 
These constraints are applied to various existing scenarios. In the case of left- or 
right-handed sneutrinos, very large electroweak gaugino masses, or very small mixing 
angles are required. 



1 Introduction 



In the cosmological Standard Model, the present day baryon density fig and the dark matter 
(DM) relic density Qdm are of similar order of magnitude: Qdm ~ 4.7 x Og. Typically one 
considers DM in the form of WIMPs (Weakly Interacting Massive Particles), which are in 
thermal equilibrium in the very early universe. Their stability can be ensured by a conserved 
quantum number of a discrete Z^v symmetry, such as -R-parity in supersymmetric theories. 
Once the temperature falls below the WIMP mass, WIMP self annihilation processes and 
the expansion of the universe reduce the WIMP density until the self annihilation rate 
falls below the expansion rate of the universe, from where on the comoving WIMP density 
remains constant. Assuming that the present day WIMP relic density coincides with Qdm 
and that its mass is in the typical range of the order 10 — 100 GeV for cold DM and for 
new particles in models beyond the Standard Model (BSM), leads to a necessary WIMP 
self annihilation cross section of the order of weak interaction cross sections. 

However, this origin of Qdm would be completely disconnected from Qb, which is con- 
ventionally attributed to a baryon asymmetry t]b originating from CP violating process 
in the early universe. Concrete calculations of Qdm show that it is quite sensitive to the 
WIMP self annihilation cross section and the WIMP mass, and varies over many orders of 
magnitude as function of the unknown parameters in BSM models attempting to explain 
the present DM relic density. Then the similar sizes of Qdm and Qb result from a numerical 
coincidence. 

Asymmetric DM (ADM) (see [IH5] for some early discussions, and [B] for a review) is 
an attempt to explain the proximity of QLdm and Qb- The particles X forming the DM are 
assumed to be distinct from their antiparticles X, and to carry a certain quantum number. 
The corresponding charge density of the universe is assumed to be related to baryon number 
through equilibrium processes in the early hot universe, such that the asymmetries r]B and 
t]x become related. If the X — X annihilation rate u X x is sufficiently large, the resulting 
X relic density will be determined exclusively by the asymmetry t]x- Then one obtains 
Qdm — ^ b (where M p is the proton mass, and M x the mass of the X particles), 

which gives the correct order if t]x ~ f]B- 

In supersymmetric scenarios, the discrete Z x symmetry responsible for the stabil- 
ity of the particles X is equal to or related to i?-parity. For instance, sneutrinos (left- 
handed or right-handed, mixtures thereof or mixtures with singlets) have been proposed as 
ADM [714T2]. Singlet extensions of the MSSM have been considered in [THHH]. Recently, 
higgsinos have been suggested in [15], since they possess a conserved quantum number be- 
fore the electroweak phase transition. In all these cases, the dark matter asymmetry can 
be related to the baryon and/or lepton asymmetry through sphaleron, gauge and Yukawa 
induced process (see [16] and references therein). As long as the baryon/lepton asymme- 
try is generated at temperatures above the freeze out of the processes which transfer it to 
the ADM, this mechanism is independent from the precise CP and baryon/lepton number 
violating origin of the baryon asymmetry. 

Of course, in the case of sizeable ADM couplings to Higgs or Z bosons, the direct 
detection rate must be studied and must not exceed present bounds [17]. (See [18] for a 
discussion of the potential conflict between a sufficiently large X — X annihilation rate and 
a too large direct detection rate.) A too large direct detection rate can be avoided through 
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mass splittings leading to inelastic scattering [121113] • Alternatively, one may consider that 
the X particles - after their density has frozen out to a value related to the baryon density 
- decay later into other essentially inert particles with very small couplings, which have 
not been in thermal equilibrium. Subsequently we will leave aside the problem of direct 
detection rates. 

However, typically the discrete Z/v symmetry responsible for the stability of the particles 
X does not forbid X — X self annihilation processes, through the same couplings which 
transfer the baryon or lepton asymmetry to the X particles. Once X — X self annihilation 
processes are allowed, these processes can wash out the asymmetry rjx- 

This conclusion is too naive: A rough condition for the absence of a wash-out is to 
require that, at temperatures of the order of the DM mass, the rate of these processes is 
below the Hubble expansion rate. Subsequently we study this phenomenon quantitatively 
in the form of the coupled set of Boltzmann equations for the symmetric and antisymmetric 
dark matter number densities. We find that the upper bounds on axx are extremely strong 
if one wishes to obtain a final X number density which is dominated by its asymmetry rjx 
such that Qdm is related to Qb as described above. 

Interestingly it turns out that, under typical assumptions as a large X — X annihilation 
rate, a fairly model independent upper bound on axx (depending in a simple way on Mx 
and the s-wave or p-wave nature of the X — X annihilation process) can be derived from the 
only condition that the X asymmetry is not reduced by a large amount. The determination 
of this upper bound is the main result of this paper. Subsequently we apply it to sneutrino 
and higgsino ADM scenarios, and to the singlet extension of the MSSM proposed in [H] 
(in an approach similar to, but slightly different from |19j). 

Our approach is based on the Boltzmann equations for the X and X number densities in 
the presence of an X — X asymmetry. Boltzmann equations in the presence of asymmetries 
have been considered previously e.g. in [T9H22] where, however, the X — X annihilation 
rate axx wa s assumed to vanish. We find that even a small X — X annihilation rate axx 
can have a strong (negative) impact on the resulting X — X asymmetry. 

In the next Section we establish the Boltzmann equation for the X — X asymmetry 
for non- vanishing axx, and clarify the assumptions allowing for its model- independent 
integration. Our main results are upper bounds on axx depending on the tolerated dilution 
of the initial X — X asymmetry. In Section 3 we study the consequences of this result for 
sneutrino and higgsino ADM scenarios. In Section 4 we consider the model with a term 
X 2 LH/A in the superpotential introduced in [2] where the expression for axx is different 
from the previous scenarios. A summary and conclusions are given in Section 5. 

2 Boltzmann equations for asymmetric dark matter 

The common mass of the DM particles X and its antiparticle X is denoted by m. Their 
equilibrium number densities are assumed to differ by a chemical potential \i. Assum- 
ing Maxwell-Boltzmann statistics (justified for temperatures T < m/3), the equilibrium 
number densities are given by 

rq = l^ gm 2 K 2 (m/T)e^ T , nf = ^- gm 2 K 2 (m/T)e-^ T . (1) 
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g denotes the number of internal degrees of freedom of X, and K 2 the modified Bessel 
function of the second kind. We assume that the X — X asymmetry has been generated 
during periods before the one considered here, through processes (such as sphaleron pro- 
cesses) which have frozen out. Since we assume that the X — X asymmetry is related to 
the baryon asymmetry, \ijT is very small: fi/T < 10~ 9 . 

It is convenient to introduce number densities per comoving volume, Y = n/s, with s 
the entropy density: 

s = he Sf (T) 2 -^T\ (2) 

where 

h eff (T)= E^(§) 3 + ^ E *(§Y ( 3 ) 

bosons fermions 

is a sum over the effective massless degrees of freedom. Furthermore one translates the 
time dependence of the number densities into a temperature dependence, using 

for the Hubble parameter in the radiation dominated epoch and 

P = 9e ff (T)^T* (5) 

for the (relativistic) matter density. g e ff is defined as in ([3]) with cubic powers replaced 
by quartic powers. Assuming a two-body final state and using the principle of detailed 
balance [23J, the Boltzmann equations for Yx and Y x as functions of x = m/T become 



a\s i — 1/2 

dYx / 7r gj m 



dx V 45G x 2 

I 1/2 

aY x / 7r g* m 



[(*xxv) {Y x - Y? 2 ) + (a xx v) {Y X Y X - Y?Yg)] , (6) 



[(a XX v) (r| - Y? 2 ) + (a xx v) {Y X Y X - Y^Yf)] (7) 



dx V 45G x 2 
where the effective number of degrees of freedom g* is given by 

Here we have assumed that the self annihilation cross sections satisfy axx = a xxi an d 
(crv) is the thermal average of the cross section times velocity given by 



(av) 



8m 4 TK 2 (m/T) 



POD 

/ a(s)(s-4m 2 ) v / ^K 1 ( v ^/T)ds. (9) 

J 4m 2 



In most cases (av) can be expanded in powers of the relative velocity of the incoming 
particles. Then, the thermal average is approximated by an expansion in powers of 1/x: 

(av) ~a + bx~ l + O (x~ 2 ) . (10) 
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It is useful to consider the difference and the sum of Y X and Y x defined by A and Z, 
respectively: 

A = Y X -Y X , Z = Y X + Y X . (11) 



Then the Boltzmann equation for A becomes 

j A I 1/2 

aA / 7T gj m 



dx V 45G x 2 
where, neglecting terms of 0(fi/T) 2 , 



c xx v) (ZA - Z"A") (12) 



Z- = A-=tz-. (13) 

27T 4 h e ff{m/x) 1 

For o xx <C o xx we can neglect a xx in the Boltzmann equation for Z, which becomes 

V2„ 



dZ it gj m 1 



dx V45G ^ te.) 5 (^-A'-Z-» + A-»). (14) 

Previously these Boltzmann equations were investigated in [191422] under the assumption 
that A remains constant, i.e. that the right hand side of ( |12p vanishes. If a xx is large 
enough, the freeze-out temperature is low, and Z ~ Z eq to a very good approximation 
over a long period, and finally Z t ^oo ~ A up to corrections studied in [T91I22] . This is the 
desired result leading to a DM relic density determined by A which, in turn, is supposed 
to be related to the baryon asymmetry. 

During the period where Z ~ Z eq , (I12D simplifies to 



which can be integrated with the usual initial condition Ai n ~ A eq for T ~ m or x ~ 1, 
and a given expression for m(a xx v). Note that it is Ai n which is assumed to be related 
to the baryon asymmetry. As in the usual case one finds that A freezes out at a freeze-out 
temperature Tf = m/xf, and = At-xx, ~ A eq (x^). However, since eq. (1151) is linear in 
A, the ratio R = A^j Ai n is independent from Ai n and hence independent from fi/T. The 
ADM paradigm requires that R is not too small; otherwise Aoo is sensitive to (cr xx v) as in 
usual DM scenarios, and Qdm ~ &b remains a numerical coincidence. 

The dominant dependence on the parameters of the model originates from the combina- 
tion m(a xx v) in (ITS"]) , an additional weak dependence on m arises from the effective number 
of degrees of freedom in h e ff(m/x) in (1131) and in g„. For these we use the parametrization 
given in [21]. Now, R can be computed for any given expression for m(a xx v). (We solved 
the coupled set fl6H7|) of Boltzmann equations, but for a xx < 10~ 5 a xx we have Z ~ Z eq for 
the relevant range of x, and the result for R is independent from a xx as it is obvious from 
(fT5"j) .) Assuming eq. ffTUl) . we will consider the two cases (a) (<J XX v) ~ a (s-wave annihila- 
tion) and (b) (<J XX v) ~b/x (p-wave annihilation; in practice one may find a combination 
of both). The results for R as function of log (ma) and log (mb) are shown in Fig. [lj 
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log (ma) or log(mb) 



Figure 1: The ratio R = A^j A in of the final asymmetry over the initial asymmetry as 
a function of log (ma) or log (ma), with ma or mb in GeV. The solid lines represents the 
values of R if {(Jxxv) = a, and the dashed lines the values of R if (axxv) = b/x. The color 
corresponds to the mass of the DM particle, red for m = 10 GeV, green for m = 100 GeV 
and blue for m = 1 TeV. 

We see that within the cases (a) or (b) the dependence on m beyond the one in m(axxv) 
is negligibly small, and we can deduce upper bounds on m(axxv) as function of the tolerated 
reduction R of the asymmetry: 

R > 0.5: ma < 5 x 10" 17 GeV" 1 (case (a)), mb < 1 x 10" 16 GeV" 1 (case (b)) , 

i? > 0.1: ma < 1 x 10" 15 GeV" 1 (case (a)), mb < 5 x 10" 15 GeV" 1 (case (b)) . (16) 

Clearly, if we require only a moderate reduction of the asymmetry A, the freeze-out 
temperature = m/x^ must not be far below m, or x^ must not be too large. (Here 
we define as the temperature where the expansion rate of the universe becomes larger 
than the rate of self-annihilation, i.e. H(x4) = (n^ — which is solved 

numerically.) In Fig. [2] we show xf as function of log (ma) and log (mb) and, indeed, xf is 
well below 5 if the first of the conditions ffl6|) is satisfied. 

The requirement of a not too large value of Xj and the fact that ((Txxv)(xf) ~ e x f 
explains why we obtain (<Jxx v ) *C (°'xx v )i the desired value of Xf in case of X — X 
annihilation is rather of (9(20). 

In the next Sections we study the implications of the upper bounds on <Jxx for various 
models for supersymmetric ADM. 
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Figure 2: The freeze-out point function of log (ma) (solid line) and log (mb) 

(dotted line). 

3 Sneutrino and higgsino ADM 

Left-handed sneutrinos or mixtures of left- and right-handed sneutrinos (or singlets) have 
been proposed as ADM in [THT2"]. Clearly, left-handed sneutrinos ul can self-annihilate 
through processes of the form vl + — > + by the exchange of electroweak gauginos 
in the t-channel. Electroweak gauginos are the binos with mass M 1; and winos with mass 
M 2 . The corresponding expression for (<J^ L ^ L v) can be obtained from [25], and reads in the 
limit Mi, M 2 > m 



where g<i is the SU(2) gauge coupling and W the weak mixing angle. If one assumes 
universal gaugino masses at the GUT scale, M\ and M 2 are related by M\ ~ M 2 /2 and, 
with tan 2 9 W m 0.3, (JXTJ) simplifies to 



From Fig. [TJ the first term leads to stronger constraints, and applying the conservative 
bound R > 0.1 (case (a)) from (fl6|) leads to 



which excludes gaugino masses of the order of the electroweak scale. (Even for such large 
gaugino masses the sneutrino-antisneutrino annihilation rate would remain large due to 
processes with slepton exchange in the t-channel.) 




(17) 




(18) 




(19) 
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If the ADM X consists in a mixture of left-handed sneutrinos Pl and right-handed 
sneutrinos or other electroweak singlets, the result depends on the Vl component of X or 
the mixing angle sin 5 where X = sin 5 + . . . : 

g\ sin 4 S ( ^ 3^1 



'""^V'-^M' < 20 > 



Now the same argument leads to 

Mo 

sin 2 S < 3.3 x 10" 6 x - . (21) 

Vm ■ 100 GeV 

Hence, for M 2 ~ m ~ 100 GeV, sin 5 has to be very small independently from constraints 
from direct DM detection experiments - or, for mixing angles sin 5 ~ 1, one is lead back 
to (H2J. 

Higgsinos h u , hd with masses m ~ 200 — 1000 GeV as ADM were proposed recently in 
|15j . Before the electroweak phase transition (where the Higgs vevs develop), higgsinos can 
be considered as mass eigenstates. Through sphaleron, gauge and Yukawa interactions at 
high temperature, a higgsino asymmetry A~ h ^ ~ A~ hd proportional to the baryon asymmetry 
is generated [151 EE]- If the higgsinos are the LSPs (lightest supersymmetric particles) 
and other sparticles are sufficiently heavy, this asymmetry could survive until today |15] 
provided the h u h u (and h d h d ) self annihilation rates are sufficiently small. 

However, higgsinos have the same couplings to electroweak gauginos as sneutrinos, and 
can again self-annihilate into H u , H d (to be considered as eigenstates before the electroweak 
phase transition) through t-channel exchange of binos and winos. This time the scattering 
process is a p-wave process (case (b) in (TTB]) ) and we obtain from [2B] (again for Mi, M 2 ^> 
m) 

3gf / tan 2 9 W 1 



— + . (22) 

Applying again the conservative bound R > 0.1 (case (b)) from (TL6]) leads to the same 
constraint on M 2 as in (|T9|) . which is close to the estimated bound on gaugino masses given 
in [15] from the same argument. 



4 The AW - XXHL/A model 

Left-handed sneutrinos and higgsinos as ADM tend to violate bounds on direct DM de- 
tection cross sections (unless mass splittings are introduced leading to inelastic scatter- 
ing [T2"|IT5]). Moreover, as we have seen in the previous section, constraints from sufficiently 
small self annihilation cross sections are very strong. These constraints would be allevi- 
ated, if the asymmetry is transferred to lighter essentially inert particles. A simple model 
of that kind has been proposed in [H], where a gauge singlet scalar superfield X and a 
superpotential 

AW = jXXH u Li (23) 

are introduced. (Here H u denotes a Higgs superfield, and Li any left-handed lepton su- 
perfield.) This nonrenormalizable interaction can originate from integrating out heavy 
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vector-like sterile neutrinos or electroweak doublets [2] with mass ~ A, typically > 1 TeV. 
Another singlet chiral superfield X should be introduced to allow for a supersymmetric 
Dirac mass term Mx for the fermionic components ^x^xi e -§- y i a a NMSSM-like singlet 
S with (S) ^ and a coupling X'SXX in the superpotential. 

ip x , ipx would carry lepton number ±1, respectively. The superpotential (|23|) breaks the 
usual R parity, but preserves a Z 4 symmetry which allows for the decay of the usual LSP 
into ipx^x |14j . At high temperature where the processes induced by ( 12 3 p are in equilibrium 
together with sphaleron, gauge and Yukawa interactions, these imply an asymmetry Ax of 
~ 35% of the baryon asymmetry [T^[T9]. the precise value depending on whether top quarks 
and squarks are still in equilibrium when the interactions from ( 123]) decouple. Assuming a 
sufficiently rapid ipx — i>x annihilation rate and Mx ~ 11 — 13 GeV, the relic density is 
then automatically of the correct order. 

After electroweak symmetry breaking, the superpotential ( 123]) gives rise to an interaction 
of the form 

jMxti ■ (24) 

(Subsequently we omit the neutrino/sneutrino index i.) The sneutrino v does not have to 
be the LSP; the LSP can be the lightest neutralino % Then an on-shell sneutrino v would 
decay via the usual vertex gvxy - where g is of the order of electroweak gauge couplings, if 
X is dominantly bino-like - into x pl us a neutrino v. At energies below the sneutrino mass 
m^, integrating out the sneutrino leads to an effective four Fermi interaction 

; "-fefeF • (25) 



mf-A 



At energies above M^, (f25|) allows for the scattering process ipxipx —> X v - However, 
assuming M% > 2Mx, the lightest neutralino % is n °t stable. Since x is a Majorana Fermion, 
([25]) leads to its decay into ipxipxv and t/jxi^xv with corresponding branching ratios of 50%. 
The latter case leads to the scattering process 

^x^x 4>x4>xvv ■ (26) 

As in the case of sneutrinos and higgsinos, this ADM self annihilation process can have 
potentially disastrous consequences for the remaining asymmetry. In [T3], the rate for 
this process has been estimated by integrating out both the sneutrino v and the lightest 
neutralino x with the result that, for ~ ~ 100 GeV and A > 1 TeV, it would go 
out of equilibrium (drop below the Hubble expansion rate) for decoupling temperatures T D 
somewhat above Mx, in which case the asymmetry would hardly be washed out. 

A more detailed analysis of the ADM self annihilation processes has been performed 
in [19]. There it was pointed out that, for ~ M^ ~ 100 GeV and A > 1 TeV, the 
dominant ADM self annihilation process is real x production through the interaction (f25]) . 
and the corresponding cross section was given. 

Subsequently the authors in [TS] estimated Tjj by equating the ADM self annihilation 
rate with the Hubble expansion rate. They tolerated a considerable wash-out of the asym- 
metry and/or a Boltzmann suppression for To < Mx, and studied the necessary relations 
between the final values for Aoo and defined in (fTTJ) (or r m = (Z^— Aoc) / (Zoo+Aoo)), the 
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DM mass Mx and the lightest neutralino mass My. required for a desired DM relic density. 
Clearly, in most of this parameter space (after a considerable wash-out of the asymmetry 
and/or Boltzmann suppression), the desired DM relic density is no longer simply related 
to the baryon asymmetry in contrast to ADM paradigm. 

Here we ask the question under which conditions this does not happen, i.e. under which 
conditions the initial asymmetry Ax determines essentially the DM relic density. As before 
we assume that the ipx — ipx annihilation rate is sufficiently large, such that we can assume 
Z ~ Z e « in (PP leading to <^E>. 

If the threshold for the process (1261) would be s > 4M X , we could apply our previous 
formulas. However, for the (dominant) annihilation process via real x (+ neutrino) produc- 
tion, the threshold is s > M| > 4M X . As a consequence the thermal average of the cross 
section times velocity in (Tl5l) depends in a more complicated way on Mx, M^ and notably 
on the temperature or x = M x /T such that the expansion ( TTUT) is no longer applicable. 
First, we use the cross section <Jxx(s) for the process ipx?Px — > XP from [19] : 

k 2 MU ( s V 

'"M^iGsr 1 ) (27) 

where k = -j^r- Then, using with M| as lower threshold of the integral (without 
expanding in M x /M^ as in [IS]), we obtain 

, s ^ / Mji V m- 4Ml)K 2 (^x) + ^KaCgfe) 

{(TXXV) = 16^ [m^ ) ^Kf(5 ' (28) 

For Mx <^ M%, (o~xxv) is proportional to Mj^ e~ M */ T like the collision term evaluated 
in [19]. The Boltzmann suppression ~ q~ m x/ t is an obvious consequence of the threshold 
s > M| required for real x production. Subsequently we integrate the Boltzmann equation 
( |T5l) numerically employing (125]) for (axxv), which allow us to study R = A^jA^ as before. 
Now, however, R depends in a more complicated way on the parameters k, My, and Mx of 
the model. On the other hand, the initial asymmetry is quite well known in this class of 
models, A in ~ 0.35 B (where B is the baryon asymmetry), and finally we must obtain 

Qdm = A^Mx ^ 4J (2g) 
Brtip 

which determines Mx in terms of or R: 

M x ~ (12.5 GeV)/i? . (30) 

Hence, once the correct ADM relic density is imposed, the free parameters of the model are 
k = j^i, M^ and Mx or R. In order to clarify the correlations between these parameters, 
we show R in the range R = 1 . . . 0.1 (Mx = 12.5 . . . 125 GeV) as function of M% for various 
values of k < 10~ 5 GeV -2 in Fig. |3j (Note that all curves continue horizontally along R = 1 
beyond their upper end.) 

The shape of these curves can be understood as follows: Near the upper end, R ~ 1, 
M x ~ 12.5 GeV and practically no wash-out takes place by construction, since (<Txx v ) is 
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Figure 3: R as function of for different values of k: red: k = 10 9 GeV 2 , green: k = 
10~ 8 GeV" 2 , orange: k = 10~ 7 GeV~ 2 , magenta: k = 10~ 6 GeV" 2 , blue: k = 10~ 5 GeV~ 2 . 
The mass Mx of the dark matter, related to R as in fl30|) . is indicated on the right vertical 
axis. 

too small. For fixed k, (crxx v ) increases for decreasing My. due to {a xx v) ~ e~ M ^ T , and 
below some critical value of M^, (cr X x v ) is large enough such that R starts to decrease. 

If a sizeable wash-out takes place (R < 0.5), it stops when the annihilation rate falls 
below the Hubble expansion rate, i.e. for a given value of (<Jxxv). For fixed k, this implies a 
certain fixed value for e~ Mi//TD or M^/Tp. Using xd = Mx/Tp and fl30|) one easily derives 
R ~ 1/(xdM^), which explains the decrease of R for increasing in this range. 

Fig. |3] allows to identify the regions in parameter space in which R is not too small 
(R > 0.5), i.e. where flnM follows naturally from as in the ADM paradigm. If the 
coefficient k relatively is large, k ~ 10~ 5 GeV" 2 or ma ~ 100 GeV and A ~ 1 TeV, one 
needs > 300 GeV. If x is lighter, ~ 100 GeV, one needs a very small value of 
k < 10~ 9 GeV~ 2 , hence correspondingly large values of ra^ and/or A. 

5 Conclusions 

If the dark matter asymmetry is related to the baryon (or lepton) asymmetry, some cou- 
plings must necessarily relate these sectors. The same couplings can lead to dark matter 
self annihilation processes, which can wash out the corresponding asymmetry A. A rough 
condition for the absence of a wash-out is to require that, at temperatures of the order of 
the DM mass, the rate of these processes is below the Hubble expansion rate. In Section 2 
we have studied the corresponding set of Boltzmann equations quantitatively (assuming a 
two-body final state). Requiring a modest wash-out (Aoo/Ai n > 0.1), the upper bounds on 
mi&xxv) are very strong, and can be deduced from Fig. [T]if (axxv) ~ a or (<r X xv) ~ b/x. 
If the ADM consists in sparticles with couplings to electroweak gauginos (left-handed 
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sneutrinos or higgsinos), it follows that the electroweak gauginos must be extremely heavy 
such that supersymmetry does not solve the hierarchy problem. If the ADM consists in 
particles like right-handed sneutrinos which mix weakly with left-handed sneutrinos, the 
mixing angle must be very small, see eq. (12 ip . 

In different models for ADM, the dominant ADM self annihilation process may be 
kinematically possible only for s above a threshold larger than (2Mx) 2 , in which case 
it becomes Boltzmann suppressed and eq. ffTUl) is no longer valid. An example is the 
popular AW ~ XXHL/A model of [Hj. Also in this case, the numerical integration of 
the Boltzmann equation for the asymmetry allowed us to specify the range of parameters 
where the wash-out of the asymmetry remains modest. 

In the past, the constraints following from the absence of a wash-out of the asymmetry 
have sometimes been neglected or underestimated (notably in the case of sneutrinos); we 
hope that the present work helps to clarify the relevance of ADM self annihilation processes 
and the resulting conditions on corresponding models. 
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